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Abstract. Measurements of the velocity of interstellar He° inside of the heliosphere have been conducted over the past forty 
years. These historical data suggest that the ecliptic longitude of the interstellar flow has increased at a rate of ~ 0.18° per 
year. Possible astronomical explanations for these short-term variations in the interstellar gas entering the heliosphere are 
presented. 
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INTRODUCTION 

When Leverett Davis, Jr., first proposed that the solar 
wind carves out a cavity in the galactic magnetic field 
[1], he noted that the path taken by galactic cosmic rays 
(GCRs) to the inner heliosphere depends on the direc- 
tions between the solar and interstellar magnetic fields 
(ISMF). Although the angle between the ISMF direction 
and the vector motion of interstellar gas flowing through 
the heliosphere was recognized as important for deter- 
mining the termination shock distance [2], the magneti- 
cally distorted heliosphere was confirmed only after Voy- 
ager 2 crossed the solar wind termination shock in 2008 
[3]. Interstellar gas and dust have now been measured 
in the inner heliosphere by many spacecraft [4], and the 
direction of the local ISMF has been independently de- 
termined by the band of energetic neutral atoms (ENAs) 
discovered by the Interstellar Boundary Explorer (IBEX) 
spacecraft [5]. The modulation of galactic cosmic rays 
(GCRs) by the heliosphere has been known since Davis 
first proposed the existence of the heliosphere, but only 
recently have asymmetries in the GeV-TeV GCRs been 
directly attributed to the modulation or acceleration of 
GCRs in the plasma heliosphere [6, 7]. 

Our heliosphere is 'blowing in the wind' of the sur- 
rounding Local Interstellar Cloud (LIC). Like all winds, 
the interstellar medium (ISM) around the heliosphere 
changes with time over both short and long timescales, 
and these variations cause changes in the flux of galactic 
cosmic rays at 1 AU that may imprint on the radionuclide 
record on Earth [8]. 

During the forty years over which the wind of inter- 
stellar He° through the heliosphere has been measured, 
the Sun has moved 150 AU through space, and 300 AU 
with respect to the LIC. Someday we might encounter 
one of the tiny AU-sized low column density interstel- 
lar "clouds", N(H°)< 10 18 - 10 19 cm" 2 , found through- 
out the ISM [9, 10]). The nearest of such clouds that has 



been identified is in the constellation of Leo, and is ~ 20 
pc away. Small scale structure of the interstellar Na° line 
indicates spatial variations in the ISM properties over 
scales of 50 AU, and cold dense tiny clouds are seen as 
well (20 K, 3000 cm~ 3 , 200 AU [11]). 

In this review, our galactic environment is briefly re- 
viewed and evidence for short-term variations in the ISM 
feeding He into the heliosphere is presented, together 
with the implications of these variations for the local in- 
terstellar environment. 



THE INTERSTELLAR MAGNETIC 
FIELD NEAR THE HELIOSPHERE 

Convergence of the predictions for properties of the outer 
heliosphere boundary conditions determined from MHD 
heliosphere models (see other papers in this volume), 
and from photoionization models of the LIC [4, 12], 
show that the heliosphere is located in a partially ion- 
ized magnetoionic medium. The discovery in 2009 of 
a 'Ribbon' of ENAs by IBEX provided an unpredicted 
and unprecedented measure of the direction of the inter- 
stellar magnetic field (ISMF) that shapes the heliosphere. 
Comparisons between 1 keV ENA data and 3D MHD he- 
liosphere models indicated that the Ribbon traces view- 
points where the ISMF draping over the heliosphere is 
perpendicular to the radial sightline [13]. The pressure 
of ENAs in the heliosheath indicates a strength for the 
ISMF shaping the heliosphere of ~ 3 /J.G [14], which is 
comparable to the expected LIC field strength if thermal 
and gas pressures are comparable [12, 15]. 

Although details of the Ribbon formation are not yet 
understood, models suggest the Ribbon originates 10- 
100 AU upstream of the heliopause [13, 16, 17, 18]. In 
the models, the locus of Ribbon points moves towards the 
equator of the ISMF as either the field strength becomes 
stronger [16, 17], or the formation region moves further 



from the heliopause [19]. 

We have undertaken a program to find the ISMF di- 
rection in the ISM near the Sun, employing measure- 
ments of the polarization of starlight that results from 
paths through clouds containing magnetically aligned 
dust grains [20]. The E- vector of the polarization is par- 
allel to the ISMF direction in the diffuse ISM (e.g. [21]). 
We developed a method to determine the ISMF direc- 
tion from an ensemble of measurements of the polar- 
ization position angles (e.g. the angle between the po- 
larization direction and the north pole). We have found 
that the nearby, < 20 pc, ISMF is aligned with the axis 
through e,b= 47° ± 20°, 25° ± 20° (optical polarization 
data do not provide the field polarity). This direction is 
32° ± 30° from the center of the Ribbon arc, located at 
^=33 ±4°,55 ±4°. 

The ISMF traced by the polarization data appears 
to have an ordered component that extends to within 
8 pc of the Sun. The ordered component of this very 
local field rotates slowly at a rate of ~ 0.25° per parsec. 
Superimposed on this field is a turbulent component of 
~ 23° [20]. Since the Ribbon arc itself can be centered 
up to 16° away from the ISMF direction outside of 
the heliosphere according to the models [16], and the 
field rotates slowly, in my view we can consider the 
ISMF directions obtained from the polarization data to 
be consistent with the direction of the ISMF that shapes 
the heliosphere. 

Figure 1 shows the direction of best-fitting ISMF di- 
rection (solid lines) against the distribution of interstel- 
lar dust within 100 pc. The large white arc-like structure 
shows the interstellar dust associated with Loop I, which 
is a giant magnetic bubble that appears to extend to the 
solar location (see below). The very local ISMF may be 
part of the Loop I magnetic field (see below). 

HELIOSPHERE AND ENVIRONMENT 

The relative Sun-interstellar velocity drives interstellar 
gas and larger dust grains into the inner heliosphere, 
while the interstellar magnetic field, ionized gas, and tiny 
charged grains are deflected around the heliopause. As 
is well known, neutral interstellar gas is ionized in the 
heliosphere by charge-exchange (CEX) with solar wind 
plasma, photoionization, and/or electron impact ioniza- 
tion, depending on the species. Charge exchange be- 
tween the solar wind and interstellar H° creates ENAs 
that are mapped by IBEX, revealing the time-variable 
heliosphere plasmas. Ionized interstellar atoms are con- 
verted outwards with the solar wind as a population 
of evolving pickup ions (PUI) that seed the formation 
of anomalous cosmic rays (ACR). All of the abundant 
atoms that have significant neutral fractions in the in- 
terstellar cloud around the heliosphere, e.g. H, He, N, 
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FIGURE 1. Directions of optical polarizations are shown 
projected against the distribution of interstellar dust within 
100 pc, as traced by color excess E(B-V) measurements. The 
polarizations trace the ISMF direction. The lines show the 
best-fitting local ISMF direction from starlight polarized in 
the ISM. The white arc-like region shows the distribution of 
interstellar dust that traces the boundary regions of Loop I 
within 100 pc. The dark regions in the third and fourth galactic 
quadrants correspond to the deficit of local ISM associated 
with the Local Bubble, which also corresponds to the regions 
with the brightest local UV and EUV radiation fluxes because 
of low interstellar opacities. The data are plotted in an aitoff 
projection, centered on the galactic center and with longitude 
increasing to the left. See [20] for more details 



O, Ne, and Ar, have been either directly detected or ob- 
served in the form of PUIs or ACRs [4] . When combined, 
these data show that the LIC is a warm, low density, par- 
tially ionized interstellar cloud, with the isotopic com- 
position of the Sun, and that arrives at the Sun from a 
direction within 15° of the galactic center [22, 4]. 

Reasonable values for the parameters of the LIC at the 
heliosphere boundaries have been reviewed in Frisch et 
al. [4]. Densities of H°, protons, and He° are, respec- 
tively, 0.19 cm- 3 , 0.07 cm" 3 , and 0.015 ±0.003 cm" 3 . 
The He density is based on Ulysses measurements near 
5 AU where particle trajectories have not yet been dis- 
torted by solar gravity into the downwind focusing cone 
[23]. Measurements of interstellar He at 1 AU by IBEX- 
LO show that ISM flows through the heliosphere at a 
velocity of 23.2 ±0.3 km s and a temperature of 
6300 ± 390 K, with the flow directed towards galactic co- 
ordinates of i,b= 185.25 o ±0.24°,-12.03 o ±0.51° (see 
Table 1 for ecliptic coordinates [24]). 

The penetration depth of neutral interstellar atoms 
into the heliosphere differs between elements, depend- 
ing mainly on the CEX cross-section. Helium and neon 
are on hyperbolic trajectories, becoming ionized within 
1 AU of the Sun, and form density enhancements down- 
wind of the Sun that are known as 'focusing tails' 
[23, 25]. In contrast, the density distribution of neutral 
interstellar hydrogen is strongly modified by CEX in the 



inner and outer heliosheath regions and inner and outer 
heliosheath regions and the heliosphere, and by radiation 
pressure closer to the Sun. The complex H° source func- 
tion, with contributions from primary and secondary par- 
ticle populations [26], temporal variations in the solar La 
flux that affect trajectories and backscattered emission 
[27], and ~ 50% filtration, challenges efforts to extract 
an accurate LIC velocity vector from the H° backscat- 
tered data. Thus the most reliable measure of the LIC 
velocity is obtained from data on interstellar He° inside 
of the heliosphere. 
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FIGURE 2. Directions of the flow of interstellar He from 
direct detection of He atoms (after 1990), and the fluorescence 
of the EUV 584A line (20th century, see Table 1). Although 
the earlier results have been discounted [28] since they did not 
agree with the results of the Ulysses GAS detector [29, 23], 
the most recent precise measurement by IBEX [24] indicate 
that there may be a real temporal variation in the direction of 
the ISM flowing into the heliosphere. If so, the results plotted 
in this figure suggest that the ecliptic longitude of the wind 
direction varies at approximately S X ~ 0.18° per year. The 
width of the horizontal bars for each point show the earliest 
and latest times over which the data were collected. The large 
blue points are the IBEX and Ulysses in situ values. 



BLOWING IN THE WIND 

Historical helium data 

The velocity vector of the flow of interstellar He 
through the heliosphere has been determined many times 
over the past 40 years. The directions of the He flow 
have been either determined from the weak resonant 
scattering of solar 584A emission by interstellar He in 
the heliosphere, or from in situ He measurements. The 
historical variations of the He° flow direction is plotted in 
Figure 2. The blue points in Figure 2 represent directions 
determined from in situ measurements of He by Ulysses 
[29] and IBEX [24]. The He velocity, beginning and 
end times of the measurements, and downwind ecliptic 



longitude of the velocity are listed in Table 1, where most 
directions have been precessed to J2000 coordinates. 

Interstellar helium enters the heliosphere with minimal 
filtration, and trajectories are determined by the position, 
velocity, and temperature of the particles, solar gravity, 
and ionization due to photoionization and electron im- 
pact ionization within 1 AU [30]. The details of the mod- 
els used to obtain the He vectors vary between the dif- 
ferent sets of data, particularly for the older data where 
the models are less developed and solar flux data were 
not reliable. 

These data provide a basis for evaluating temporal 
variations in the interstellar boundary conditions of the 
heliosphere. Figure 2 shows that the ecliptic longitude 
of the He velocity vector has steadily increased with 
time over the past forty years (Table 1). A linear fit 
to the data suggests an increase of ~ 0.18° per year 
in the ecliptic longitude of the He flow. Over the same 
interval the ISMF rotation would be only l/100th of an 
arcminute based on the polarization data. At this rate, 
the direction of the LIC flow would change by 90° over 
500 years. Depending on the behavior of the ISMF, the 
heliosphere configuration should also vary dramatically 
over timescales of several centurys.. 

Assuming that the time-dependence of the He direc- 
tion is real, I will now look at possible interstellar expla- 
nations for this variation. 



Why the change ? 

It has been argued that the older data are inaccurate 
because of instrumental uncertainties or poorly known 
ionization corrections. In an effort to make the older 
Prognoz 6 He° data [35] consistent with the Ulysses 
GAS results [29], Lallement et al. [28] found that the 
two sets of results could be brought into agreement if the 
instrumental noise was decoupled from the pre-launch 
laboratory calibration and added as a variable into the 
analysis. However this type of argument can not explain 
the even steeper temporal variation between the wind 
longitudes found by IBEX versus Ulysses (Figure 2). 

One might speculate that the slow increase in ecliptic 
longitude of the velocity is due to variations in solar con- 
ditions over the past 40 years. At a given spatial position, 
the He° atom may have arrived directly from the ISM, 
or may have arrived indirectly through close passage by 
the Sun on the downwind side due to gravitational focus- 
ing [39]. An apparent shift in the direction could occur 
if there were a systematic variation in the mix of direct 
and indirect atoms that are detected, perhaps due to vari- 
ations in the ionization rates of the indirect atoms trav- 
eling close to the Sun where losses from electron impact 
ionization and photoionization dominate. Such an expla- 
nation would not explain the difference between Ulysses 



TABLE 1. Interstellar gas velocity obtained from interstellar He data 
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Ref. 


1972.8/1973.6 


5-20 


73.1 ±3, 


-8.34±3 


2000 


Weller & Meyer [31] 


1974.1/1974.1 


[22 ±3]** 


73.4±5, 


-5.4±5 


2000 


Ajello et al. [32] 


1975.6/1976.1 


9-27 


73. ±3, - 


-7.34±3 


2000 


Weller & Meier[33] 


1976.9/1977.1 


22-28 


73.6±3, 


-4.37 ±3 


2000 


Weller & Meier [34] 


1977.8/1978.1 


27 ±3 


75.2±3, - 


- 6 ± 3 xxx 


2000 


Dalaudier et al. [35] 


1990.9/2002.7 


26.3 ±0.4 


75.4±0.5, 


-5.1±0.2 


2000 


Witte [29] 


1992.6/1993.6 


26.4 ±1.5 


76.0±0.4, 


-5.4±0.6 




Flynn et al. [36] 


1992.6/2000.9 


24.5 ±2.0 


74.7±0.5, 


-5.7±0.5 




Vallerga et al. [37] 


2009.2/2010.2 


23.2±0.4 


79.00 ±0.47, 


-4.98±0.21 


2000 


McComas et al. [24] 



* Column 1 shows the first and last year of data acquisition. 

1 J2000 coordinates were either obtained by precessing B1950 celestial coordinates, or originally provided as 
J2000 coordinates. Other values were taken directly from the references. 
** The Copernicus velocity was assumed from Adams & Frisch [38]. 



and IBEX values, since some Ulysses data were acquired 
near 5 AU before the flow characteristics become dis- 
torted by solar gravity. Sunspot activity levels vary be- 
tween the different sets of measurements, but no system- 
atic trend is associated with the data listed in the Table. 

A more interesting possibility, that I prefer, is that the 
direction of the ISM flow through the heliosphere is time- 
dependent. 



Cloud edges 

Cloud edges are interesting properties of the ISM that 
are difficult to study, and the Sun is near a cloud edge. 
It is now well known that an absorption component at 
the LIC velocity has not been seen towards stars in the 
upwind direction. The most compelling example is 36 
Oph, which is a star 6 pc away and 9° from the helio- 
sphere nose direction. Towards this star, the velocities of 
the interstellar Fe+, Mg+, and D° lines are 28.3 ±0.5 
km s _1 [40] in contrast to the ~ 23.2 km s _1 that is pre- 
dicted from the He° LIC velocity. This indicates that for 
an average LIC density of 0. 19 cm~ 3 , the Sun is less than 
20,000 AU (0. 1 pc) from the LIC boundary. 

The properties of the LIC cloud edge depend on the in- 
teraction of the LIC and the adjacent ISM. Within 10 pc, 
only 25%-40% of space is filled with observable ISM. 
The LIC is part of a flow of ISM (the cluster of local 
interstellar clouds, CLIC) with an upwind direction di- 
rected towards the center of the Loop I evolved super- 
bubble, more specifically the S 1 shell part of Loop I and 
the flow is decelerating ([4, 41], and see below). A possi- 
bility for the edge region of the LIC in the upwind direc- 
tion is that it is influenced by cloud collisions, since the 
flow of nearby ISM is not a rigid body but rather is com- 
pressed towards the Sun from both the extreme upwind 
and extreme downwind directions [42, 43]. For example, 



Schwadron et al. [6] explain the GCR asymmetries re- 
lated to the heliosphere in terms of a magnetic field dis- 
turbance caused by the collision of the LIC and another 
cloud in the downwind direction. A second possibility is 
that the LIC is embedded in the hot plasma of the Lo- 
cal Bubble, and a hot conductive interface or turbulent 
mixing layer is present at the LIC edge. 

The one piece of information that we have is that the 
change in the flow direction is such that the flow is mov- 
ing away from the direction of the ISMF that shapes the 
heliosphere, which is towards A,j3= 221°, 39° accord- 
ing to the center of the IBEX Ribbon arc [5, 44]. Note 
that the latitude variations are smaller and less well con- 
strained than the longitude variations. This would seem 
to indicate that the ISM feeding the He into the he- 
liosphere has become less ionized over time, since in- 
creased ionization would couple the flow more tightly to 
the ISMF. The decreasing-ionization scenario is consis- 
tent with a local gradient in the interstellar radiation field 
that is defined by the lower opacity and higher radiation 
fluxes found in the galactic interval £~ 190° - 360° (Fig- 
ure 1, [20]. 



Turbulence 

Based on the trend for the ecliptic longitude of the LIC 
velocity vector to increase with time, it should be con- 
sidered a real possibility that the Sun has sampled turbu- 
lence in the LIC velocity over the past several decades. 

The average values for the temperature and turbu- 
lence of the LIC are 7, 500 ± 1,300 K and 1.62 ±0.75 
km s _1 , based on absorption lines toward 19 stars be- 
hind the LIC according to Redfield et al. [45]. Nonther- 
mal turbulence is calculated by assuming that the full- 
width-half-max of the absorption line is proportional to 

(thermal + turbulence) where ^thermal IS the Doppler 



spread in thermal velocities appropriate for the cloud 
temperature, and Vturbuience represents non-thermal con- 
tributions to the velocity over the length of the sightline. 
Uncertainties in this type of analysis can incorporate gra- 
dients in the velocity field, or even blended absorption 
components, into the turbulence parameter depending on 
spectral resolution and data quality (e.g. [46]). 

If the average line width of the LIC is attributed to 
the ISM entering the heliosphere, but with a temperature 
<~ 6, 300 K as found from the He data, then the turbulent 
velocity would be 4.7 km s _1 . There is no information 
about the scale size of the turbulence in the LIC, so 
according to these arguments, a plausible case can be 
made that this nominal 4.6 km s -1 turbulent component 
samples turbulence in the ISM between the Sun and the 
upwind edge of the LIC. In such a case, a systematic 
variation in the flow direction is possible. 

Filaments in the magnetoionized medium 

A gradient in the magnetic field strength between the 
heliosphere and the edge of the LIC is also possible. Tur- 
bulence in the galactic magnetic field manifests itself in 
the magnetoionic ISM over large spatial scales through 
structures that appear in maps of the polarized radio con- 
tinuum, but which are not visible in the radio contin- 
uum intensity maps. Gaensler et al. [47]. have shown that 
these structures form a web of filaments that do not ap- 
pear to be supersonic. Most interesting for the local ISM 
is that the spatial gradient of the radio continuum polar- 
ization is such that the polarization changes more rapidly 
for directions that are perpendicular to the filamentary 
elongations. If the LIC is elongated along the ISMF di- 
rection as given by the IBEX Ribbon, the gradient in the 
ISMF would make an angle of ~ 42° with respect to the 
He° flow vector. The gradient in the ISMF would influ- 
ence the heliosphere and secondary He feature discov- 
ered by IBEX [48], and possibly the mix of primary and 
indirect He° atoms that are detected. 



Loop I superbubble 

My last possible explanation for the variable ISM di- 
rection inside of the heliosphere is a hand-waving argu- 
ment, based on the dynamic structure of the ISM near the 
Sun, which seems to be related to the Loop I superbub- 
ble. Three epochs of star formation in the nearby Sco- 
Cen Association formed Loop I (e.g. [4]). Comparisons 
between optical and synchrotron polarization data show 
that the ISMF is parallel to the elongated filaments that 
form the observed structure of Loop I (e.g. [49]). The 
brightest part of Loop I in the 1.8 MHz radio continuum 
is a region known as the North Polar Spur, where the den- 



sity of the ISM increases continuously from the Sun out 
to - 100 pc [50, 51, 20]. The local ordered ISMF dis- 
cussed above is directed towards the North Polar Spur. 

Three sets of data suggest that the Sun is near or in a 
fragment of the shell of the Loop I superbubble: (1) The 
Loop I model of Wolleben [52, 53] identifies two radio 
continuum shells (called "SI" and "S2") at different 
distances and locations. The Sun is located in the rim of 
the S 1 shell according to his model, a result that agrees 
with prior estimates of the Loop I configuration [54] . (2) 
Interstellar clouds within about 20 pc of the Sun, e.g. 
the CLIC, have a bulk velocity in the local standard of 
rest (LSR, e.g. after solar motion with respect to nearby 
cool stars is subtracted) that is within 15° of the center of 
the S 1 shell, suggesting an expanding shell configuration 
[53, 4]. (3) The local ISMF directions found from the 
center of the IBEX Ribbon arc, and from the interstellar 
polarization measurements, both have angles ~ 76° away 
from the bulk local ISM flow, and ~ 64° from the SI 
center as defined by Wolleben. 

In conclusion, there are many astronomical justifica- 
tions for the hypothesis that the very local ISM varies 
over spatial scales that are comparable to the heliosphere 
size. After all, the heliosphere itself provides a major 
disruption of the ISM and ISMF over spatial scales of 
> 10,000 AU, or equivalently 0.05 pc. A reanalysis of 
the historical He° data is needed to check my basic con- 
clusion that the direction of the flow of interstellar gas 
through the heliosphere has varied over the space age. 
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